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Abstract: We review the status of CPT violation in the neutrino sector. Apart from 
LSND, current data favors three flavors of light stable neutrinos and antineutrinos, with 
both halves of the spectrum having one smaller mass splitting and one larger mass splitting. 
Oscillation data for the smaller splitting is consistent with CPT. For the larger splitting, 
current data favor an antineutrino mass-squared splitting that is an order of magnitude 
larger than the corresponding neutrino splitting, with the corresponding mixing angle less- 
than-maximal. This CPT- violating spectrum is driven by recent results from MINOS, but 
is consistent with other experiments if we ignore LSND. We describe an analysis technique 
which, together with MINOS running optimized for muon antineutrinos, should be able 
to conclusively confirm the CPT- violating spectrum proposed here, with as little as three 
times the current data set. If confirmed, the CPT- violating neutrino mass-squared differ- 
ence would be an order of magnitude less than the current most-stringent upper bound on 
CPT violation for quarks and charged leptons. 
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1. Introduction 

All known particles are either self-conjugate under CPT or have CPT conjugate "antipar- 
ticles". In every case the antiparticle partner is observed to have the same mass as the 
corresponding particle, within experimental resolutions. These observations are consistent 
with the description of all non-gravitational particle interactions by local relativistic quan- 
tum field theory, where CPT conservation is a result of the intimate connection between 
Lorentz invariance, locality, hermiticity, and the absence of operator ordering ambigui- 
ties. For precisely this reason it is important to pursue increasingly rigorous tests of CPT 
invariance, and to extend our experimental constraints to sectors previously beyond reach. 

In this regard neutrinos are especially interesting. Neutrinos have tiny nonzero masses, 
suggesting that the neutrino mass generation mechanism has novel features and that neu- 
trinos communicate to a sector of new physics whose effects on charged leptons and quarks 
are as yet unobservable. As demonstrated in the next section, the current generation of 
neutrino oscillation experiments are sensitive to CPT- violating effects orders of magnitude 
smaller than what so far could have been detected for charged leptons or quarks. There is 
both theoretical and experimental motivation to pursue a rigorous study of CPT properties 
for neutrinos, keeping in mind that CPT violation may correlate with other exotic effects 
such as Lorentz violation or quantum decoherence. 

In this report we update [l]-[8] the experimental constraints on CPT violation for 
neutrinos, focusing on the case where other new physics effects are subdominant to a 
CPT- violating difference in neutrino/antineutrino mass spectra. As favored by the data 
we also assume three flavors of light stable neutrinos and antineutrinos, both halves of the 
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spectrum having one smaller "solar" mass difference and one larger "atmospheric" mass 
difference. For the larger splitting we show that the global data set favors an antineutrino 
mass-squared splitting that is an order of magnitude larger than the corresponding neutrino 
splitting, as well as an antineutrino mixing angle 623 that is less than maximal. This CPT- 
violating spectrum is driven by recent results from MINOS [9] , but is consistent with other 
experiments. 

We describe an analysis technique to confirm or deny the best- fit CPT- violating hy- 
pothesis with future data. We advocate and demonstrate the use of the Neyman-Pearson 
hypothesis test [10], also known as the a- (3 test, generalized from ratios of simple likeli- 
hoods to ratios of extended likelihoods with floating parameters. This method has the 
advantage, for a given likelihood ratio, of distinguishing between the test significance a, 
the probability that CPT-conserving masses and mixings are rejected even though they 
are in fact correct, and the power of the test 1 — f3, where (3 is the probability that the 
CPT- violating solution is rejected even though it is in fact correct. For a given future data 
set, one can require that the p- value of the CPT-conserving hypothesis as extracted from 
the likelihood ratio is less than some benchmark significance a chosen according to one's 
theoretical prejudice about CPT violation. 

In advance of new data we can use Monte Carlo experiments to extract the value of (3, 
thus estimating the prospects for distinguishing CPT violation in the neutrino spectrum 
if it is in fact present. We examine these prospects for the MINOS experiment. To be 
conservative, in maximizing the likelihoods we do not float parameters defining the CPT- 
violating mass spectrum, since this would tend to increase the maximum likelihood for 
the CPT- violating hypothesis even when it is wrong. We do however float experimental 
parameters related to the overall neutrino production rate and the energy spectrum; float- 
ing these parameters increases the maximum likelihoods for the incorrect hypotheses while 
leaving the maximum likelihoods for the correct hypotheses essentially unchanged, thus 
lessening the power of the Neyman-Pearson test. 

Even with this conservative approach, we demonstrate that MINOS running optimized 
for muon antineutrinos should be able to conclusively confirm the CPT- violating spectrum 
proposed here, with as little as three times the current data set. 

2. CPT violation in the neutrino sector 

2.1 CPT violation with and without Lorentz violation or other exotic new 
physics 

The discovery of parity (P) violation in fundamental interactions was a big surprise, espe- 
cially considering that P is an element of the extended Lorentz group. As we now under- 
stand, it is possible to violate P in quantum field theory without compromising invariance 
under the restricted Poincare group that includes only proper orthochronous Lorentz trans- 
formations, i.e. Lorentz transformations continuously connected to the identity. 

For CPT, the connection to Lorentz invariance is even stronger. As emphasized by 
Feynman, in a local description of quantum field theory the Lorentz invariance of off-shell 
amplitudes requires combining processes with propagation of both off-shell states and CPT 
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conjugates of those states. Going the other way, Greenberg has shown [11] that in quantum 
field theory CPT- violating mass differences on-shell inevitably lead to Lorentz-breaking 
effects off-shell, with consequences for both locality and operator-ordering in quantum 
field theory. 

Because of the initimate theoretical connection between CPT and Lorentz invariance, 
experimental searches for CPT violation are related to experimental tests of Lorentz in- 
variance. In both cases the most straightforward experimental approach is to look for 
departures from the expected relativistic on-shell dispersion relations for particles and an- 
tiparticles: 

2 - *2 2 — 2 — * 2 q / \ 

E = p + m , E =p +m , m = m , (2-1) 

where here and throughout a bar denotes a quantum number of a CPT conjugate state. 
This relationship suggests three experimentally distinct scenarios: 

• Detectable violations of Lorentz invariance in the dispersion relations for some par- 
ticles, but conserving CPT to within experimental resolutions. 

• Detectable violations of Lorentz invariance in the dispersion relations for some par- 
ticles, accompanied also by detectable violations of CPT. 

• Detectable violations of CPT in the dispersion relations for some particles, but con- 
serving Lorentz invariance to within experimental resolutions. 

The first two scenarios are motivated by the possibility of a spontaneous breaking of vacuum 
Lorentz invariance, perhaps related to new Planckian physics such as space-time foam, su- 
perstrings, or extra dimensions [12]-[22]. The third scenario is motivated by the possibility 
of non-local physics whose primary on-shell effect may be CPT violation [23]. 

A further complication is that exotic new physics such as quantum decoherence [12]-[14] 
or extra dimensions [21] may lead to baseline-dependent effects on neutrino oscillations with 
additional CPT- violating features not captured by deviations from the expected dispersion 
relations. Of course matter effects, though predicted by the Standard Model, are also an 
example of baseline-dependent effects on neutrino oscillations with CPT- violating features. 

For neutrino oscillation experiments there are thus effectively three kinds of tests of 
CPT: 

1. Searches for Lorentz-violating effects in concert with CPT violation. The current 
best limits on this case are from the MINOS experiment [24]; we will not elaborate 
further on this scenario. 

2. Searches for CPT- violating differences between neutrino and antineutrino mass spec- 
tra. This is the main subject of our report. 

3. Searches for inconsistencies in oscillation results that could signal baseline-dependent 
new physics with possible ramifications for CPT. Prospective limits are discussed in 
[16, 19] and [21]. 
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Figure 1: The CPT- violating neutrino spectrum proposed in [2] as an explanation of LSND. 

In the last case there is an important connection between CP and CPT. Even when 
CPT is conserved, CP violation in neutrino mixing allows the possibility of differences 
between neutrino and antineutrino oscillation probabilities in neutrino appearance experi- 
ments: 

P(y a n) + P{y a u h ) • (2.2) 

However, as shown in [4], CP violation without CPT violation cannot produce a neutrino- 
antineutrino discrepancy in disappearance experiments: 

P{Va Vfi) ¥= P(Pa Vi) ■ (2.3) 

A corollary of these results is that a neutrino-antineutrino oscillation discrepancy arising 
from CP violation without CPT violation requires at least two relevant mass splittings 
contributing to the oscillation, as occurs e.g. in some (3+2) sterile neutrino models [25]. 

2.2 Comparing limits on CPT violation 

Assuming that the source of CPT violation is a mass asymmetry in the dispersion relations 
2.1, the relevant figure of merit in comparing different experimental limits on CPT- violation 
is the mass-squared difference between a particle and its CPT conjugate. 

For quarks the most stringent experimental limit [26] is from neutral kaons, whose 
mass-squared difference is constrained to be less than 0.5 eV 2 , or ~ 0.1 eV 2 if we attribute 
the CPT asymmetry to the constituent strange quarks. 
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Figure 2: The 3+2 CPT-conserving but potentially CP-violating neutrino spectrum proposed 
in [25] attempting to reconcile LSND, MiniBooNE, and other short-baseline oscillation results. 



For charged leptons, the most stringent constraint [26] is from the upper limit on the 
electron-positron mass difference; this corresponds to an upper bound on the mass-squared 
difference of approximately 2 x 10 4 eV 2 . 

The CPT- violating best fit reported here corresponds to a difference of mass-squared 
differences of only 0.02 eV 2 . This means that for neutrinos the current generation of 
oscillation experiments have sensitivities to potential CPT- violating effects orders of mag- 
nitude smaller than the above limits. Note that Bahcall et al. reached the same conclusion 
applying different figures of merit [27] (see also [28, 29]). 

Thus, contrary to what is sometimes implied in the literature, it is plausible that 
CPT- violating mass differences would be detected first in the neutrino sector, even if such 
effects have comparable magnitude in the quark and charged lepton sectors. Furthermore, 
as noted already in the introduction, since neutrinos appear to gain mass through a novel 
mechanism, it is also plausible that CPT- violating mass differences are much larger for 
neutrinos compared to the other sectors. 
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3. CPT-violating neutrino mass spectra 

The main constraints on the mass differences and mixings of neutrinos come from the 
neutrino oscillation experiments [30]- [44] summarized in Table 1. 
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Table 1: Summary of current and past neutrino oscillation experiments. The first column shows 
the principle oscillations that the experiment could in principle observe; the second column indicates 
whether this constitutes an appearance or disappearance experiment. The third column indicates 
the primary sensitivity either to "solar" mass splittings, "atmospheric" (atm), or "short baseline" 
(SBL). 

Because we are interested in the possibility of CPT violation, we will consider the 
masses and mixings of the neutrino mass matrix as completely independent of the the 
masses and mixings of the antineutrino mass matrix, and consider the experimental con- 
straints on each matrix separately. Because of the flavor sensitivity of the SNO results, 
the active neutrino composition of the "solar" neutrino oscillation is well-constrained. The 
Super-Kamiokande data also have some flavor sensitivity; technically this measures the 
sum of the "atmospheric" neutrino and antineutrino oscillations, but in practice is mostly 
constraining for the neutrinos, which dominate over the antineutrinos as cosmic ray secon- 
daries in the relevant energy range. The Super-K atmospheric neutrino data are bolstered 
by accelerator-based experiments K2K and MINOS, which report muon neutrino disap- 
pearance consistent with the atmospheric mass splitting and large mixing. The net result 
[45] is that the active neutrino masses and mixings are required to closely resemble the 
left half of the spectrum shown in Figure 1, modulo the possibility of inverting the solar- 
atmospheric hierarchy. The main question on the neutrino side is whether there are small 
admixtures of one or more sterile neutrinos in the three light mass eigenstates, but as yet 
there is no evidence for such mixings. 

On the antineutrino side, the situation is less clear. KamLAND has reported an 
electron antineutrino disappearance signal consistent with an antineutrino mass splitting 
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and mixings equivalent to the "solar" counterpart on the neutrino side. LSND reported 
a Vn —* v e appearance signal consistent with an antineutrino mass-squared splitting ~ 1 
eV 2 . MINOS has reported preliminary muon antineutrino disappearance results, consistent 
with an antineutrino mass-squared splitting that is roughly the geometric mean of the 
KamLAND and LSND favored splittings. 

Thus, even allowing for CPT violation, oscillations between three active antineu- 
trino species cannot reconcile KamLAND, MINOS and LSND simultaneously. The CPT- 
violating spectrum shown in Figure 1, proposed in [2] to accommodate solar, atmospheric, 
and LSND splittings with only three active flavors, was conclusively excluded by Kam- 
LAND [6]. 

Without resorting to new baseline-dependent exotic physics, this leaves two possibili- 
ties: 

• Case (i) The LSND results are incorrect. 

• Case (ii) The LSND results are correct, but the corresponding short-baseline (SBL) 
oscillation involves mixing with one or more species of sterile neutrinos. 

In the second case one may question whether it is even necessary to resort to a CPT- 
violating mass spectrum, since the addition of sterile neutrinos adds new parameters that 
potentially loosen up the experimental constraints. Several recent analyses [25],[46]-[48] 
have looked at this question in detail, using global fits that (in the case of [48]) include the 
latest MiniBooNE data. The conclusion is that CPT-conserving sterile neutrino scenarios, 
even allowing for the possibility of large CP violation in the case of two or more sterile 
species, cannot avoid at least a 3 a discrepancy among different experimental data sets, with 
the largest tension between the SBL appearance experiments and the SBL disappearance 
experiments. 

Thus Case (ii) requires either that we disregard some oscillation results other then 
LSND, or that we again resort to a CPT-violating spectrum. Thus for example one could 
develop CPT-violating versions of the 3+2 spectrum shown in Figure 2. We will not pursue 
this possibility further here, since it is already under investigation elsewhere [49]. 

The remainder of this paper is devoted to Case (i): we disregard the LSND signal, 
and explore to what extent a CPT-violating neutrino spectrum is allowed, or even favored, 
by the remaining global data set. For simplicity we will assume that the situation is not 
further complicated by sterile neutrinos or baseline-dependent exotic physics, though of 
course both are possible. 

With these assumptions there is qualitatively only one CPT-violating neutrino mass 
spectrum candidate, shown in Figure 3 and first discussed in [6]. To be more precise there 
are four candidate spectra, since we can invert the hierarchy on either the neutrino or 
antineutrino side independently, but existing data is insensitive to these choices, with the 
exception of the neutrino observations from supernova SN1987A [1]. 

As discussed above, the neutrino side of the spectrum in Figure 3 is completely con- 
strained by data. On the antineutrino side, the smaller "solar" mass splitting is necessary 
to accommodate the KamLAND signal; a CPT-violating variation of this splitting and 
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the related mixings is still allowed at about the 5% level. The larger mass splitting has 
to accommodate the antineutrino disappearance signals from MINOS and Super-K, the 
null appearance results from KARMEN and MiniBooNE, as well as the null disappearance 
results from other SBL oscillation experiments. 
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Figure 3: The best-fit CPT-violating neutrino spectrum obtained from our analysis. 



4. Constraining the antineutrino spectrum 

To make detailed contact with the experimental results we first introduce the neutrino 
survival and transition probabilities given by 



P(v a -> vp) = 5 af 3-4 ^2 U ai UpiU a jUf3j sin 2 
i>j=i 



for neutrinos and 



P(u a -> up) = 5 a p-A ^2 U ai U piU aj U pj sin z 

%>j=i 



4E 



(4.1) 



AE 



(4.2) 



for antineutrinos. The matrix U = {U a i} (U = {U a i}) describes the weak interaction 
neutrino (antineutrino) states, i> a , in terms of the neutrino (antineutrino) mass eigenstates, 
i/j. That is, 



yiUqjVi and 



(4.3) 
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where we have ignored the possible CP phases. The matrices can be parametrized as 
follows: 




C12C13 S12C13 
U = I -S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 I (4.4) 
S12S23 - C12C23S13 -C12S23 - S12C23S13 

and similarly for U. In Eq. (4.1) L denotes the distance between the neutrino source and 
the detector, and E is the lab energy of the neutrino. 

We use the notation Am s 2 D | ar = Am 2 2 , Am 2 tm = A7n 2 3 to denote the smaller and 
larger mass-squared splittings on the neutrino side, and Am s 2 okr = Arnf 2 ; Am 2 tm = Am 2 3 
for the antineutrinos. 

SBL reactor experiments give important constraints on the antineutrino spectrum. 
Their results indicate [30, 31, 32] that electron antineutrinos produced in reactors remain 
electron antineutrinos on short baselines. Because of the short baselines we can ignore 
the smallest ( "solar" ) antineutrino mass difference and average the other two; the survival 
probability can be expressed as 

P(V e ^V e ) = l-2U%(l-U%). (4.5) 

Thus, even for rather large antineutrino mass differences, the survival probability will be 
close to one if is either almost one or almost zero. Physically this means that we can 
choose between having almost all the antielectron flavor in the heavy state (which really 
means the furthest away state since we can invert the spectrum) or alternatively leave this 
state with almost no antielectron flavor. The first possibility was depicted in the Figure 1, 
while the second is realized in Figure 3. 

MINOS and Super-Kamiokande constrain both the larger antineutrino mass-squared 
difference Am 2 tm and the antineutrino mixing angle 623- KamLAND constrains mostly the 
smaller antineutrino mass-squared difference Am s 2 olar . 

We have performed a x 2 fit °f the antineutrino spectrum (assuming three active flavors 
only) using the data from MINOS, Super-Kamiokande, KamLAND, and CHOOZ. The 
best-fit result is shown in Figure 3. The CPT- violating features are encapsulated in: 

Am 2 tm = 0.02 eV 2 , sin §23 = 0.407 , (4.6) 

compared to the global-fit neutrino spectrum values 

Am 2 tm = 0.0025 eV 2 , sin 23 = 0.707 . (4.7) 

This CPT violation is driven by the MINOS results; indeed our best-fit values for 
Am 2 tm and sin 2 2#23 are close to those reported by MINOS fitting their data alone. 

The overall quality of our fit is good, with a x 2 P er degree of freedom of 0.98. As 
seen in Figure 4, the x 2 deviation as a function of the single variable Am 2 tm has a clearly 
defined minimum. We note however that this is only the case when #23 is allowed to float 
in the fit; if 623 were fixed to maximal mixing, then the chi-squared distribution in Am 2 tm 
would be rather flat. 
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Figure 4: The relative \ 2 deviation of the CPT-violating fit to data from MINOS, Super- 
Kamiokande, KamLAND, and CHOOZ, as a function of the single parameter A?7i^ tm (in eV 2 ), 
with sin 623 = 0.407. The lower (red) curve excludes the MINOS muon antineutrino data, while the 
upper (blue) curve includes it. 

5. Discussion and future prospects 

The MINOS muon antineutrino disappearance results should be regarded as preliminary. 
They are from data runs with the target optimized for neutrinos, introducing more com- 
plicated systematics for the antineutrinos, and poorer statistics (42 events observed at the 
far detector). This situation will improve dramatically with results from MINOS running 
optimized for antineutrinos, scheduled to begin soon. 

Our fit shows that large, order-of-magnitude CPT violation in the neutrino sector 
is still a viable possibility. Making the further assumptions that the CPT violation is 
(approximately) baseline-independent and does not have a strong dependence on sterile 
mixing, a unique CPT-violating mass and mixing pattern is selected, up to the four-fold 
ambiguity of inverting the neutrino and/or antineutrino hierarchies. 

The most timely question is whether better data in the near future from MINOS 
could provide compelling evidence for neutrino CPT violation. To address this question, 
we have performed toy muon antineutrino disappearance experiments, using the survival 
probability obtained either from a CPT-conserving spectrum or from our best-fit CPT- 
violating spectrum. We use the reconstructed muon antineutrino energy spectrum reported 
by MINOS, but to add some realism we allow a one-parameter distortion of the energy 
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spectrum, and float this parameter in the fit. We also float iVo, the mean expected number 
of neutrinos detected in the MINOS far detector in the absence of oscillations; while this 
number is estimated in the experiment, it is subject to significant systematic uncertainty. 
We also float N osc , the actual (but unmeasured) number of neutrinos in each experiment 
that would have been detected had they not oscillated to a different neutrino flavor. 
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Figure 5: The log likelihood ratio distribution for 300,000 toy experiments simulating 200 nominal 
muon antineutrinos per experiment on the MINOS baseline. The right-hand histogram uses the 
CPT-conserving hypothesis to generate the toy results, while the left-hand histogram uses the 
CPT- violating hypothesis. In both cases the ratio is from the maximum likelihood computed with 
the CPT-conserving pdfs over the maximum likelihood computed with the CPT- violating pdfs. 
The histograms are normalized to unit probability. 

Each toy experiment is equivalent to MINOS running with 200 nominal muon antineu- 
trino events expected in the far detector in the absence of oscillations. This is approximately 
a factor of 3 increase over the current data. For each of 300,000 toy experiments based 
on each mass spectrum, we compute the maximum likelihood {i.e. we maximize the likeli- 
hood with respect to the floated parameters) for both CPT-conserving and CPT-violating 
hypotheses; then we plot the normalized distribution of events versus the logarithm of the 
ratio of the likelihoods. The definition of the likelihood and the details of the analysis are 
presented in the appendix. The result is shown in Figure 5. 

This plot allows a Neyman- Pearson test of the CPT-conserving versus CPT-violating 
hypotheses. We choose a cut a on the log ratio of the likelihoods for the case that the 
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toys are based on the CPT-conserving spectrum; the significance a corresponds to the 
probability that the CPT-conserving hypothesis is rejected even though it is true. Clearly 
we should choose a small value for a, since we have a strong prior bias that CPT is 
conserved. Having thus fixed a we can extract (3, the probability that the CPT-conserving 
hypothesis is accepted even though the CPT- violating spectrum is the correct one. Then 
1 — (3 is the measure of the power of this hypothesis test. 

The results are very encouraging: even with a chosen as small as 6 x 10~ 7 , correspond- 
ing to a Gaussian significance of 5a, we find 1 — j3 very close to unity. This indicates a 
nearly 100% chance that the CPT-violating spectrum discussed in this paper would be 
correctly chosen by the hypothesis test if it is in fact true. 

Figure 6 shows the oscillation probabilities of both the CPT-conserving and CPT- 
violating hypotheses plotted as a function of the muon antineutrino energy. The current 
MINOS data points (binned in energy) are superimposed. From this figure it is clear 
that the CPT-violating hypothesis makes clear energy-dependent predictions about what 
should be observed in future MINOS running: 

• The lowest energy bin will rise. 

• The dip apparent to the eye around 10 GeV will remain. 

This figure also explains why our prediction for the power of Neyman-Pearson test 
with 200 nominal MINOS events is so encouraging, even though we have very conservatively 
floated the total number of neutrinos in the likelihood fits. The discrimination of the CPT- 
conserving and CPT-violating hypotheses comes both from a significant difference in the 
overall survival probabilities and from the dramatic difference in the energy dependence of 
the survival probabilities. Even allowing for a rather large systematics, as we have done 
here, data generated from one hypothesis is almost never as well-described by the incorrect 
hypothesis, for experiments with at least 200 events. 
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A. Likelihood methods and neutrino oscillations 

Consider a typical neutrino oscillation disappearance experiment, in which rii neutrinos of 
a particular flavor are observed in a far detector in some number of energy bins labelled 
by i. Let N s = ^T raj be the total number of neutrinos observed, and for simplicity ignore 
the possibility of fakes or neutrinos from background sources. 

Using observations in a near detector or some other method, one computes the mean 
number of neutrinos Nq that one would have expected to observe in the far detector in the 
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Figure 6: Comparison of the binned MINOS muon antineutrino disappearance data with the 
oscillation curves obtained from the CPT-conserving hypothesis (blue) and the CPT-violating 
hypothesis (red). 

absence of neutrino oscillations. For simplicity assume that Nq is the mean of a Poisson 
distribution, although one could also handle more complicated distributions. The experi- 
ment will also calculate the energy distribution of neutrinos expected at the far detector 
in the absence of oscillations. Both of the aforementioned distributions are hypotheses, 
perhaps with floating parameters representing uncertainties, but in both cases assume that 
the hypothetical distributions are independent of the particular neutrino oscillation model 
being tested. 

Now suppose one has a hypothesis for the correct neutrino oscillation model. This 
amounts to specifying the pdf p s (E), the probability that a neutrino of energy E does 
not oscillate to a different flavor. Convolving these pdfs with the energy distribution one 
obtains p l s , the probability that a neutrino is in the ith energy bin and did not oscillate. 
Letting p l ° tal = ^2 i p i s , the probability that a given neutrino does oscillate to a different 
flavor is just p osc = 1 — p*°* a '. 

The appropriate binned extended likelihood function given all these assumptions is 
given by 

_ e -"°(jVo)^" (N s + N osc )l „ „ 

(N s + N osc )l N s !N osc ! Posc i-\ {Ps) ' (AJ) 

i=i 

where N n * r denotes the total number of neutrinos that oscillated to a different flavor. 



- 13 - 



In most applications of the extended likelihood formula, the analog of the total number 
of events, here N s + N osc , is known, while the mean expected number Nq is obtained from 
the fit by maximizing the likelihood. In the case at hand an estimate of Nq is already given, 
iV s is measured, and N osc is unknown. Thus one would like to obtain N osc by maximizing 
the likelihood. From the explicit dependence shown in A.l, it is easy to see that the 
likelihood is maximized by solving 

M^osc + 1) = log[iVo Pose ] , (A-2) 

where tpQ is the digamma function. 

For N osc 2, an excellent approximation to the solution of A. 2 is given by 

N osc = N oPosc - - . (A.3) 

Substituting back into A.l, we obtain an expression for the likelihood function already 
maximized with respect to the floating value of N osc : 

N s \T[N oPosc + l] Posc l^ 8 * V ' 

This likelihood can then be further maximized with respect to other floating parameters. 

In a real experiment the mean expected number Nq is estimated from other data with 
some error. If we take this error to be Gaussian, we can include this distribution in the 
definition of the extended likelihood, and maximize the likelihood with respect to both Nq 
and N osc . The likelihood function becomes 

exp (- (JV o~5 0) ) p -N 0(N \N s +N 0Posc -l/2 ^ 
L = ^ ^_J_^_ [M ' p^oPo.c-1/2 TT( ^ (A . 5) 

where the parameters A^o and aj\r are supposed to be fixed by, e.g., extrapolating from 
near detector data. In our fits we have used Nq = 200 and <tjv = 20. We then maximize 
the likelihoods allowing Nq to take any value such that N osc is nonnegative. This increases 
the maximized likelihood of the wrong hypothesis while having very little effect on the 
maximized likelihood for the correct hypothesis; thus floating the value of Nq decreases the 
log likelihood ratio of the correct hypothesis over the wrong hypothesis. 

In our analysis we introduced a parameter c to represent uncertainty in the normalized 
energy distribution of neutrinos reaching the MINOS far detector: 

e -aE £b 

P(E) = q _ (6+1) F[b - y , a = 0.193489; b = 1.43356 + c, (A.6) 

where the numerical constants come from a fit to the MINOS spectrum. By allowing c to 
vary from —0.2 to 0.2 in the fit independently for each neutrino oscillation hypothesis, we 
introduce a variability in the energy spectrum as illustrated in Figure 7. This increases 
the maximized likelihood of the wrong hypothesis while having very little effect on the 
maximized likelihood for the correct hypothesis; thus allowing a distorted energy spectrum 
decreases the log likelihood ratio of the correct hypothesis over the wrong hypothesis. 
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Figure 7: The one-parameter variability in the normalized muon antineutrino energy spectrum 
floated in our maximum likelihood fits. The central (blue) curve is our best fit to the binned MINOS 
spectrum. The horizontal axis is energy in GeV. 
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